Fine tuning of phase qubit parameters for optimization of fast single-pulse readout 
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We analyze a two-level quantum system, describing the phase qubit, during a single-pulse readout 
process by a numerical solution of the time-dependent Schroedinger equation. It has been demon- 
strated that the readout error has a minimum for certain values of the system's basic parameters. In 
particular, the optimization of the qubit capacitance and the readout pulse shape leads to significant 
reduction of the readout error. It is shown that in an ideal case the fidelity can be increased to 
almost 97% for 2 ns pulse duration and to 96% for 1 ns pulse duration. 
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In recent years superconducting Josephson junction 
circuits have attracted a considerable interest as promis- 
ing devices for quantum computations To speed- 
up the readout of a qubit state and improve its fidelity, 
the single-pulse readout technique has been suggested [Sj- 
(Io| . According to this technique, a measurement of a 
flux-biased phase qubit, described by a shallow potential 
well with two energy levels |0) and |1) (see Fig. [1]), is 
performed by subjecting the qubit to a pulse driving. It 
leads to lowering the barrier between the qubit " left" and 
"right" potential wells, so the system will tunnel from 
state |1) with a probability close to one, while state |0) 
remains intact. However, for small readout times the 
achieved fidelity was rather low, of order 70-80 %, which 
was explained by different sources of decoherence. Re- 
cently, it has been demonstrated that the coherent Rabi 
oscillations remain nearly unaffected by thermal fluctua- 
tions up to temperatures of IK [ll[ (i.e., until the energy 
of thermal fluctuations kT becomes comparable with the 
energy level spacing Tiu of the qubit), so without degrad- 
ing the already achieved coherence times, phase qubits 
can be operated at temperatures much higher than those 
reported so far. This may signal that relatively large 
readout errors of practical devices [(| can be attributed 
to non-optimal readout of the qubits rather than quan- 
tum and thermal fluctuations. In 12J it has been demon- 
strated that the readout error has a minimum as function 
of both amplitude and duration of the readout pulse, as 
well as the qubit capacitance. On one hand it allows to 
obtain high fidelity for fixed pulse duration by designing 
the qubit with a proper capacitance, and by changing the 
pulse amplitude. On the other hand it is possible to de- 
crease the readout time, making it much smaller than the 
qubit coherence time and suppressing the effect of ther- 
mal fluctuations [Hj]. As it has been understood 
the pulse shape has a strong influence on different error 
probabilities. Also, it is believed that the main source 
of error during the qubit readout is due to incomplete 
discrimination between the two quantum states, which 
is determined by the depth of a shallow potential well. 
However, the investigation versus external magnetic field 
(which is changing the potential well depth and can be 
varied during an experiment), and also existence of the 
optimal readout pulse shape were not studied. The main 
goal of the paper is to find the optimal values of qubit 
parameters for very fast readout of order 1 ns. 




FIG. 1: The profile of a bistable potential. Solid curve - the 
original potential, dashed curve - the potential with enlarged 
deep well to simulate the effect of damping. Inset: the con- 
sidered trapezoid-like pulses with walls, from top to bottom, 
sin 2 (87rf/£ p ), sin 2 (47rf/t p ), sin(27rf/t p ), sva 2 {2-Kt/t p ). 



We describe a phase qubit by_ 
(shown in Fig. 1, solid curve) [11 



the following potential 



V(x,t) =Ej{{x-^p{t)f/2i- 
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Here Ej = Ic<h/2e is the Josephson energy, x is the 
Josephson phase, e is the electron charge, and h is the 
Planck constant. The qubit parameters are taken the 
same as in [8j]-floT]: the critical current Iq = 1.7/iA, 
the inductance of the ring L = 0.72 nH and the capac- 
itance C = 700 fF, correspond to I = 2eIcL/h = 3.71, 
2e 2 /hC = 0.6933 x 10 9 Hz, Ej/h = I c /2e = 5.31 x 10 12 
Hz. It is convenient to introduce the "inverse capaci- 
tance" D = 2e 2 /hC x 10~ 9 Hz, and express the time in 
nanoseconds. The dimensionless external magnetic flux 
tp(t) = 2ir[ciQ + Af (t)] consists of two components: the dc 
component an , which determines the depth of the shallow 
well, and the driving readout pulse with an amplitude A 
and pulse shape f(t), the trapezoid-like function, which 
grows and drops by different laws (see the inset of Fig. [T]) . 
We note that t p is defined as the full width of the pulse at 
zero level, rather than the full width at half maximum. 

The number of the discrete energy levels can be char- 
acterized by the crude estimated value: 



(2) 



where AUi is the depth of the left well (the energy 
difference between the potential maximum and mini- 
mum) and LJi is the classical oscillation frequency near 
the left-well bottom (the plasma frequency): oji = 
y/Ej(l/£ + cos xi)/m, where xi corresponds to the left- 
well bottom and m — %/ (2D) is the effective mass. 

Let us consider the readout error N, which is the sum 
of two probabilities, P±q not to tunnel during the pulse 
action from the state |1), and Poi to tunnel from the state 
|0) (i.e. N = Pio+Pqi, while the fidelity F = 1-N). The 
investigation is performed via computer simulation of the 
Schrocdinger equation and is focused on the readout error 
N versus the pulse amplitude and the shape, as well as 
the depth of a shallow potential well. 

The Schroedinger equation for the wave function 
*S>(x,t) has the following form: 



N'> 
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hC dx 2 



*(z,t). (3) 



The Eq. ([3]) does not take into account either damping 
or noise, and in the frame of Eq. ([3]) we will only study 
the error occuring due to incomplete discrimination of 
the quantum states |0) and As it has been done 



12| . to prevent the repopulation error, we consider a 



modified potential, which differs from Eq. (fTJ) that at the 
bottom of the deep well the potential does not grow up 
and remains constant to the right really far away, see Fig. 
1, dashed curve. Therefore in the boundary conditions 
*(c, t) = *(d, t) = we take c = -3, d = 797, while the 
right- well minimum is located at x r ~ 6 and xi s» 1.4. 
We have checked that further increase of this range does 
not change the results. 

For fast single-pulse readout system we take the pulse 
duration t p — 2 ns. A special feature of the quantum 
system is the impossibility to use rectangular readout 
pulses fl3 . which in the classical case leads to minimal 



noise-induced errors [13( . In the qubit a rectangular pulse 
leads to nonadiabaticity of the tunneling event, which 
results in lifting to higher eigenstates and considerable 
increase of the probability to tunnel from |0) state and 
thus to much larger values of N than for all other pulse 
shapes. As it will be shown, to obtain the minimal error 
we should use a compromise shape, close to a meander, 
but not yet leading to the discussed effect. 

First, let us take a sine-trapezoid function which grows 
and drops for t < t p /8 and t > 7t p /8 by sin 2 (47rf/i p ). 
The readout error N versus the pulse amplitude A is 
presented in Fig. 2 for different values of the inverse ca- 
pacitance D and the shift of a potential barrier ao = 0.81. 
It is seen that N has a minimum depending on the ampli- 
tude, as well as depending on the parameter D. Decrease 
of D — s- leads to increase in the number of energy levels 
and to decrease in the distance Huji — ¥ between them 
(see Eq. @, where ~ V~D), which complicates the 
discrimination between the states |0) and |1). If huji is 
large, both ground and excited states become too close 
to the barrier top, so the tunneling from both states may 
occur even without the driving pulse, leading to large 
readout error. Therefore, since both limits of large and 
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FIG. 2: The readout error N versus the pulse amplitude A 
for the pulse with walls sin 2 (47rf/i p ) and oo = 0.81. 



small D lead to large readout errors, there must be some 
optimal value of the inverse capacitance, leading to the 
minimal readout error N(A,D). 

Note that to find the minimal N(A) for a fixed value 
of D it is not necessary to plot the whole curve N(A). 
Changing the inverse capacitance leads to a shift of the 
curves on a specific number A s , e.g. increase of D by 0.3 
leads to a decrease N m i n (A) on A to 0.0037. Thus, one 
can predict the location of N(A) minimum for different 
D 1 that specds-up cither calculations or measurements. 
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FIG. 3: The readout error N versus the pulse amplitude A 
for the pulse with walls sin 2 (Ant / t p ) with different inverse 
capacitances D and dc magnetic field ao- 

The barrier height AUi depends on the external flux 
tp(t), particularly it is determined by the dc component 
a . If the well is too shallow, the variation of the pa- 
rameter D does not seriously change the small number 
of the discrete levels Ni (Eq. and the tunneling er- 
ror is still large. If this well is deep enough, the value 
Ni remains large and discrimination between two nearby 
states |0) and |1) is complicated. The optimal curves 
with minimal readout error N m i n (A, D) for different val- 
ues of external shift ag are presented in Fig. 3. It demon- 



strates the absolute minimum N m i n (A, D, ao)=0.031 for 
a =0.77, D=l.9 Hz, .4=0.0625 (fidelity F=0.969). The 
location of this minimum corresponds to some value be- 
tween five and six energy levels inside the shallow well, 
which is larger than for ao=0.81 in [13]. The fast qubit 
state measurement requires that the maximum external 
flux ipi — 27r[ao + A] changes the potential so that the 
state |1) is very close to the barrier top. It corresponds 
to JVj slightly larger than unity. So we can facilitate the 
optimization process finding a range of A values for the 
fixed parameters ao and D. For example, for ao=0.82 and 
D=l.2 Hz the minimal error corresponds to the values of 
A from 0.0166 (for the condition N t < 1.1) to 0.0236 (for 
the condition Ni > 1). While for a =0.7 and D=2 Hz, 
Nmin is shifted to A from 0.131 to 0.139. These results 
are confirmed by numerical solutions shown in Fig. 3. 
Note that this estimate is valid not only for pulse shapes 
considered in this paper, but even for the triangular and 
linear ramp pulses. 

The second pulse we consider is the sine-trapezoid 
function which grows and drops for t < t p /4 and t > 
3i p /4 by sin 2 (27rt/t p ). Using the developed algorithm, 
the curves with minimal readout error N m i n (A, D) were 
found for different constant magnetic field components 
ao (Fig. 4, solid curves). For comparison the optimal 
curves for the pulse with walls sin 2 (47rt/t p ) and the same 
ao (with their best inverse capacitance D) are also pre- 
sented by dashed curves. Relative to the previous pulse, 
in this case the absolute minimum shifts toward smaller 
values of ao and larger A, while the tunneling error N 
increases to about 0.034. So, reducing the flat part of 
the sine-trapezoid function leads to worse results. 
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FIG. 4: The readout error N versus the pulse amplitude A 
for the pulse with walls sin 2 (2vrt/t p ) with different D and ao 
(solid curves) . The curves for the pulse with walls sin 2 (47rt /t p ) 
with the same ao and their own D are given by dashed curves 
for comparison. 

Let us increase the width of the flat part of the pulse 
to 14t p /16 value. The readout error N versus ampli- 
tude A for dc magnetic field ao=0.81 is presented in Fig. 
5 (solid curves - pulse with walls sin 2 (87rf/£ p ), dashed 
curves correspond to sin 2 (47rt/t p ) for the same parame- 



ters). As expected, if the pulse shape is more close to 
a rectangular one, the effect of exitation of |0) state is 
more noticeable: there are several local minima of N (for 
D = 1.1, N = 0.044, and for D = 2.1, N = 0.045). But 
the absolute value of the error in this case is larger than 
for the sine-trapezoid pulse with walls sin 2 (47rf/t p ). 
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FIG. 5: The readout error N versus the pulse amplitude A 
for the pulse with walls sin 2 (8nt/t p ) and ao = 0.81 (solid 
curves) . The dashed curves correspond to the pulse with walls 
sm 2 (4-Kt/t p ) with the same values of D and ao- 

Thus, we have found that for the sine-trapezoid pulse 
sin 2 (47rf/i p ) with the width t p —2 ns the readout error can 
be reduced to 7V=0.031 and this value is obtained for the 
inverse capacitance £>=1.9 Hz. Therefore, considering 
the fixed set of experimental parameters of [ij]-[l(3|, the 
reduction of C from 700 fF to the range from 450 fF 
to 150 fF must lead to the highest fidelity. It has been 
demonstrated 12j that, with decrease in duration, the 
minimum N is shifted towards larger A, but the value of 
tunneling error increases significantly. For t p =l ns, sine- 
trapezoid function with walls sin 2 (87rf/t p ), ao = 0.77, A 
= 0.057 and D = 3.4, we achieve N w 0.038, which leads 
to the increase of the readout error in about 20 % only 
in comparison with 2 ns pulse. Without the pulse shape 
and dc flux optimization for t p —2 ns N was 0.053 [l2|, 
and the difference of the errors for 1 ns and 2 ns durations 
was about 40-50 %. 

Finally, let us consider a task, where the qubit capaci- 
tance C is selected out of the optimal range. In this case 
one can set some different pulse shape to maximize the 
fidelity. For example, let us take the inverse capacitance 
as in 0-0, £"=0.6933, corresponding to (7=700 fF. 
Here the highest fidelity F=94.8 % is achieved by using 
smoother pulse shapes, such as sine-trapezoid pulses with 
walls sin 2 (27rf/£ p ) and sin(27rf/t p ), while sin 2 (Ant / t p ) 
leads to larger readout error. 

In conclusion, we have performed computer simula- 
tions of the fast single-pulse readout process of a two- 
state quantum system. It has been demonstrated that 
the minimization of the readout error can be achieved 
by variation of the depth of a shallow potential well of 
a qubit. Considering the concrete parameters of existing 



qubit designs [8| - [lOj , it is recommended to decrease the 
qubit capacitance down to 150-450 fF. Further improve- 
ment of the qubit fidelity can be reached during an ex- 
periment by adjustment of the dc magnetic field and the 
readout pulse amplitude. It is demonstrated that there 
is an optimal pulse shape minimizing the readout error, 
which for C — 450 fF is close to sine-trapezoid function 
with walls sm 2 (Airt / t p ) . For larger values of the capaci- 
tance smoother pulses can lead to the maximal fidelity. 
Finally, the performed optimizations allowed to reach al- 
most 97% fidelity for 2 ns pulses and to 96% fidelity for 1 
ns pulses. While in our model several important sources 
of decoherence are not taken into account, the obtained 
values of the fidelity seem to be the highest observed for 
so short pulses (in particular, the readout error was de- 
creased by a factor of 2 in comparison with Ref. so 
we believe that the described optimizations can help to 
improve the fidelity of real devices. 

The work was supported by RFBR (projects 09-02- 
00491 and 08-02-97033). 



* Electronic address: alp@ipm.sci-nnov.ru 



[1] Y. Makhlin, G. Schon, and A. Shnirman, Rev. Mod. 
Phys. 73, 357 (2001). 

[2] A.N. Korotkov, Quantum. Inf. Process. 8, 51 (2009). 

[3] J.M. Martinis, Quantum. Inf. Process. 8, 81 (2009). 

[4] Yu.A. Pashkin, O. Astafiev, T. Yamamoto, Y. Nakamura, 
and J.S. Tsai, Quantum. Inf. Process. 8, 55 (2009). 

[5] O. Buisson, F. Balestro, J. P. Pekola, and F.W.J. 
Hekking, Phys. Rev. Lett. 90, 238304 (2003). 

[6] K.B. Cooper, M. Steffen, R. McDermott, R.W. Sim- 
monds, S. Oh, D.A. Hite, D.P. Pappas, and J.M. Marti- 
nis, Phys. Rev. Lett. 93, 180401 (2004). 

[7] O. Astafiev, Yu.A. Pashkin, Y. Nakamura, T. Yamamoto, 
and J.S. Tsai, Phys. Rev. Lett. 93, 267007 (2004). 

[8] R. McDermott, R.W. Simmonds, M. Steffen, K.B. 
Cooper, K. Cicak, K.D. Osborn, S. Oh, D.P. Pappas, 
J.M. Martinis, Science, 307, 1299 (2005). 

[9] Q. Zhang, A.G. Kofman, J.M. Martinis, and A.N. Ko- 
rotkov, Phys. Rev. B 74, 214518 (2006). 
[10] A.G. Kofman, Q. Zhang, J. M. Martinis and A.N. Ko- 
rotkov, Phys. Rev. B 75, 014524 (2007). 
[11] J. Lisenfeld, A. Lukashenko, M. Ansmann, J.M. Martinis, 

and A.V. Ustinov, Phys. Rev. Lett. 99, 170504 (2007). 
[12] A.L. Pankratov and A. S. Gavrilov, Phys. Rev. B 81, 
052501 (2010). 

[13] A.L. Pankratov and B. Spagnolo, Phys. Rev. Lett. 93, 
177001 (2004); A.V. Gordeevaand A.L. Pankratov, Appl. 
Phys. Lett. 88, 022505 (2006). 



